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Abstract 

Using 232 fb^ 1 of data collected by the BABAR detector, the E' c + and E'® baryons are reconstructed 
through the decays: E c + — ► E+j and E c ° — ► E®j, where — ► E~^ + t^ + and E® — ► E~ir + . By 
measuring the efficiency-corrected yields in different intervals of the center-of-mass momentum, the 
production rates from B decays and from the continuum are extracted. For production from B 
decays, the branching fractions are found to be B{B — > E c + X) x B{E^ — > E~tt + 7t + ) = [1.69 ± 
0.17 (exp.) ± 0.10 (model)] x 10~ 4 and B(B -» E' C X) x B{E° -» E~tv+) = [0.67 ± 0.07 (exp.) ± 
0.03 (model)] x 10~ 4 . For production from the continuum the cross-sections are found to be 
a(e + e~ -> e'+X) x B(E+ -> ~-ir + ir + ) = 141 ± 24 (exp.) ± 19 (model) fb and a(e + e~ -> E' C X) x 
B(E® — ► S , ~7r + ) = 70 ± 11 (exp.) ± 6 (model) fb. The helicity angle distributions of E' c decays are 
studied and found to be consistent with J = ^. 
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1 INTRODUCTION 



Charmed baryons have a complex and intricate spectroscopy. The states can be classified in the 
following categories: 



State 


Quark content 


Isospin configuration 


A c 


cud 


Isosinglet 


E c 


cqq 


Isotriplet 


— c 


esq 


Isodoublet 




CSS 


Isosinglet 



where q indicates a u or a d quark. There are numerous states for each of these quark flavor 
configurations , several of which have been observed [2] . 

In this analysis, we focus on the S c states, and in particular on the lowest resonances above the 
ground state, the E' c . This is one of three E c states, listed below, which are not radially excited 
and which have zero orbital angular momentum: 



State 


Approx. mass (MeV/c 2 ) 


Light flavor wavefunction 


J p 




2470 


Antisymmetric 


i + 

2 


'c 


2580 


Symmetric 


1 + 
2 




2645 


Symmetric 


3 + 
2 



Note that the J p of the E' c and E* have not been directly measured but are assigned from the 
quark model predictions. In a recent study of the decay E® — > fl~K + the helicity angle of the 
S 1 ^ was found to be consistent with J = ^, though higher spins were not excluded. 

The mass difference between the E' c and ground state is only Am = m(E c ) — m(E c ) — 
107 MeV/c 2 . Hence, the E' c is below threshold for a strong decay via S c tt or other final states 
such as A C K or DA, and so it can only decay by photon emission, E' c — ► E c ^. 

The E c + and E c ° states were observed by CLEO in 5.0 fb -1 of data [I]. This observation has not 
yet been confirmed by another experiment. Charmed baryons can be produced in two ways at e + e~ 
B-factories: from the continuum and in decays of B mesons. The CLEO measurement was made 
with requirements that x p > 0.5-0.6 depending on the decay channel, where x p = p*/ys/4 — m%, 

and p* is the center-of-mass momentum of the S' c . These requirements suppressed combinatorial 
background, but they also removed any E' c production from .B-decays, retaining only the continuum 
production of S' c . 

Recent results from BABAR and Belle indicate that B decays to charmed baryon pairs occur at 
a high rate even when the available phase space is small El HI El E] • One possible explanation is 
that baryon formation is enhanced when the two baryons are almost at rest in their center-of-mass 
frame, and strongly suppressed when their relative motion is large (see Ref. J01 and the references 
therein). This may also explain threshold enhancements seen in several modes such as B~ — > ppK~ 
and B~ -> ppK° [TT1IT2]. If processes with low energy release are favored, the production of low- 
lying charmed baryon resonances such as S' c may be substantial. There is currently no experimental 
evidence for the production of S' c in B decays. 

2 THE BABAR DETECTOR AND DATASET 

The data used in this analysis were collected with the BABAR detector at the SLAC PEP-II asym- 
metric energy storage ring. A total of (231.9 ±3.5) fb _1 of data are used, of which 210.3 fb _1 were 
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taken on the T(4S) resonance (y/s = 10.58 GeV) and the remaining 21.6 fb 1 were taken below 
the BB threshold (yfs = 10.54 GeV). The on-resonance data contains (228.3 ±2.5) x 10 6 BB pairs. 
The BABAR detector is described elsewhere |13j . 

Simulated events with S' c decaying into the desired final states are generated for the processes 
e + e~ — > cc — > E' C X and e + e~ — > T(4S) — > BB — > S' C X, where X represents the rest of the 
event. The pythia simulation package [12], tuned to the global BABAR data, is used for the cc 
fragmentation and for B decays to S' c , and geant4 ^3] is used to simulate the detector response. 

3 ANALYSIS METHOD 

3.1 Overview 

The S' c + and E'J 3 are reconstructed through the following decays: 

where the daughter S c baryons are reconstructed as follows: 

S~ — ► Air~ 

A — ► p7T _ . 

The measured invariant mass spectra of the S£ candidates are corrected for efficiency. From the 
S' c yields in different p* intervals, the production rates from B decays and from the continuum are 
extracted. We also study the helicity angle (Oh) distribution. 



3.2 Selection and reconstruction of and 

A A candidate is reconstructed by identifying a proton with dE/dx and Cherenkov angle measure- 
ments |13j and combining it with an oppositely charged track interpreted as a tt~ , and fitting the 
tracks to a common vertex. The A candidate is then combined with a negatively charged track 
interpreted as a tt~ , and fitted to a common vertex to form a S~ candidate. For each A and S~ , 
the invariant mass is required to be within 3a of the central reconstructed value, where a is the 
fitted mass resolution (approximately 1.0 MeV/c 2 for A and 1.5 MeV/c 2 for The invariant 

mass is then constrained to the nominal value [2|. Each resulting S~ candidate is then combined 
with one or two positively charged tracks interpreted as 7r + to form a £® or candidate. No 
mass constraint is applied to the S c candidates. 

Since the S~ has a long lifetime (cr = 4.9 cm), we improve the signal-to-background ratio by 
rejecting prompt background. The displacement vector from the event primary vertex to the S~ 
decay vertex is required to be at least 2.5 mm in the plane transverse to the beam direction. In 
addition, the scalar product of the displacement vector with the S~ momentum vector is required to 
be positive, rejecting unphysical candidates. These criteria were optimized in a previous analysis [J] 
and were finalized before the m(S c 'y) spectrum was examined in data. 

The invariant mass distributions for the S c candidates satisfying these criteria are shown in 
Fig. E The fitted function (black curve) is the sum of two Gaussian functions with a common mean 
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Figure 1: Invariant mass spectra of the E c ground states in 232 fb 1 of data. The reconstructed 
candidates of (a) via — ► E~tt + tt + , and (b) E® via E® — > E~ir + are shown. 

for the signal plus a first-order polynomial for the background. The half-widths at half-maximum 
(cr) of the signal lineshapes are 7.0 MeV/c 2 and 7.5 MeV/c 2 for and E®, respectively. 

3.3 Reconstruction of S c + and S'® 

Clusters of energy in the electromagnetic calorimeter are identified. The cluster must spread over 
at least two crystals. Clusters which lie along the trajectory of a charged track in the event are 
eliminated. The remaining clusters define photon candidates. The energy of the photon candidate 
must be at least 30 MeV and the lateral moment (defined in Ref. ^5]) mu st be less than 0.8. The 
E c candidates shown in Fig. ^ are then combined with each of the photon candidates to form E' c 
candidates. The mass difference Am is then computed: 



Am 



m(E 7r + 7r + 7) — m(E tt + tt + ) for E' c + 
m(S~TT +/ j) — m(E~ir + ) for E'® 



Since Am and m(E~7r + [ir + ]) are essentially uncorrelated and the mass difference between E' c and 
E c is small, this method gives good mass resolution and avoids the need to apply a mass-constraint 
to the E c candidates. To make the spectra easier to interpret, we plot the mass difference Am plus 
a constant offset m fj- sct which is approximately equal to the nominal ground state mass; offsets 
of 2.467 GeV/c 2 and 2.471 GeV/c 2 are used for the and E® states, respectively. The spectra 
thus obtained are shown for E^ and E® in Fig. [2j For each mode, the mass spectrum is shown for 
p* > 0.0 GeV/c (upper), p* > 2.5 GeV/c (middle), and p* > 3.5 GeV/c (lower). Clear E' c signals 
are observed at a mass of approximately 2.58 GeV/c 2 . 

The photon selection criteria are quite loose, especially in comparison to the previous CLEO 
analysis which imposed a minimum photon energy of 100 MeV in addition to requirements on 
the lateral shower profile and a veto of photons from 7r° candidates. The reason for the different 
selection strategies is evident from a comparison of the energy spectra of photons from S' c produced 
in the continuum, E' c produced in B decays, and background (Fig.EJ: the CLEO study excluded E' c 
from B decays and was optimized for sensitivity to continuum production of E' c where the photon 
energy is a powerful discriminant between signal and background, whereas the selection criteria for 
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Figure 2: The S^y invariant mass spectra, shown with the following p* requirements: 0.0 GeV/c 
(a,b), 2.5 GeV/c (c,d), 3.5 GeV/c (e,f). The left column shows S+7 and the right column shows 
£^7. The shaded histograms are taken from the S c mass sidebands (5a-8a from the central value, 
where a is the S c mass resolution), and the solid points are from the S c signal region (within 3a 
of the central value). 
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Figure 3: Photon energy distributions. The dotted and dashed histograms show photons from 
simulated E' c decays where the E' c is produced in B decays or from the continuum, respectively. 
The solid histogram shows photons from E' c candidates in the data with m(E~ tt + [tt + ]) in the E c 
mass sidebands and Am in the signal region (0.09 < Am < 0.12 GeV/c 2 ). 

this analysis were chosen to retain E' c from B decays with high efficiency. For illustrative purposes, 
the mass spectra of E' c candidates with tighter selection requirements are shown in Fig. 

3.4 Contributions to the invariant mass spectra 

Four principal categories of events contribute to the E' c candidate distributions: 

1. Signal E' c — > E c ^ decays which peak in both m(E~ir + [ir + ]) and Am. This is shown in 
Fig. El (a). 

2. Combinatoric background which does not peak in either m(E~Tr + [ir + ]) or Am. This is shown 
in Fig. (b). 

3. Background where a real E c is combined with an unrelated photon candidate. This peaks in 
m(E~ir + [ir + ]) but not in Am. This is shown in Fig. 0(c). 

4. Background contribution from events where a real E' c — ► S c j decay occurs and the correct 
photon is found but the E c is partially mis-reconstructed. This is shown in Fig. 0(d). This 
category generally does not peak in m(S~Tr + [Tc + ]) but peaks in Am (provided the momentum 
of the fake E c candidate is close to the real E c momentum). 

Categories 2 and 3 do not peak in Am, so we describe them with a smooth polynomial function. 
The Am distribution of the fourth category is almost indistinguishable from the signal distribution. 

One further possible contribution to the mass spectrum was considered: feed-down from the 
decay E* — ► E c tt° where only one of the two photons produced in the 7r° decay is used, leaving the 
same E c j final state as for a E' c decay. A study of this process with a simple kinematic simulation 
indicates that it would produce a very broad, non-peaking structure in the E' c mass spectrum, and 
therefore falls into the third category {E c background not peaking in Am) already discussed. 
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Figure 4: The invariant mass spectra, requiring that p* > 2.5 GeV/c, that the photon energy 
be above 200 MeV, that the shower contain at least two crystals, that the lateral moment be less 
than 0.6, and that the shower be well-contained (shower energy within two cells of the maximum 
be at least 90% of the shower energy within one cell of the maximum). Plot (a) shows S^j and 
plot (b) shows Sj?7. The mass windows are narrower than for Fig. |2j the shaded histograms are 
taken from the S c mass sidebands (5a-7a from the central value) and the solid points are from the 
S c signal region (within 2a of the central value) . 

3.5 Fitting procedure 

The data are divided into ten p* intervals of width 0.5 GeV/c from 0.0 to 5.0 GeV/c. For each p* 
interval, the Am distributions are fitted with the combination of a signal lineshape extracted from 
the simulated signal events and a second-order polynomial function to describe the background. 
The signal lineshape is parameterized as the sum of three Gaussian functions, parameters of which 
are determined from a fit to a high-statistics sample of simulated signal events in the corresponding 
p* range. The lineshape is described in more detail in Section [A. II of Appendix EI During this fit 
all nine parameters of the triple Gaussian function are allowed to vary independently within fixed 
ranges. 5 The lineshape is then fixed for fitting the data. 
The fit to the data is performed in several steps: 

1. First, we fit just the background function to an upper mass sideband of S' c , 2625 < Am + 
"^offset < 2900 MeV/c 2 , using a binned maximum likelihood method. This provides the initial 
values of the background parameters. 

2. Next, we attempt to fit the combined background and signal functions in the mass range 2550- 
2700 MeV/c 2 . The signal mass and yield are floated, with the initial value of the signal yield 
set to zero. The background parameters are also left free; all other parameters are fixed. The 
mass is initially set to a central value extracted from a fit to the entire dataset and is allowed 
to vary within 10 MeV/c 2 around this value. The fit uses a binned maximum likelihood 
method, followed by a binned % 2 minimization with MINOS error-handling enabled |17j . 

5 Fewer signal events were generated at the extremes of the p* spectrum. In cases where fewer than 1,000 signal 
events were reconstructed, only two Gaussian functions were used. This applies to 0.0-0.5 GeV/c for , and to 
4.0-4.5 GeV/c for both and 3°. For 4.5-5.0 GeV/c no signal events were generated, so the fits of 4.0-4.5 GeV/c 
were used. 
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Figure 5: Illustrations of the signal and background contributions for 5^7 in simulated continuum 
events. Plot (a) shows category 1, correctly reconstructed signal. Plot (b) shows category 2, 
combinatoric background. Plot (c) shows category 3, background with real S c . Plot (d) shows 
category 4, where a photon from a real S' c decay is combined with an incorrectly reconstructed S c . 
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3. We then check whether the fit converges to a physical value. If the signal mass is within 
2 MeV/c 2 of the edge of the allowed range or the yield is unphysical (above 50,000 or be- 
low — 100), we reject the fit. This typically occurs when studying a region of phase space 
where the signal yield is too small to fit with a floating mass. 

4. If the first fit is rejected, we reset the parameters to the initial values described in step 2 and 
fix the mass to the central value extracted from the fit to the entire dataset. The fit is then 
repeated. 

The individual fitted spectra are shown in Appendix IE! 

To remove the category 4 background described in Section f3. 41 we perform a sideband subtrac- 
tion in m(S^TT + [TT + ]) as follows: 

• The Am distribution of events in the S c mass signal region — 3a < m(S~ir + [ir + ]) — ttiq < +3a 
is plotted, where a is the S c mass resolution and mo is the central value of the S c mass peak. 

• Similarly, the Am distribution for events in the S c mass sidebands, 5a < \m(S~7r + [ir + ]) — mo| < 
8a, is plotted. 

• The Am distributions are fitted with a signal lineshape plus a polynomial background as 
described above. The integral of the signal function gives the combined yields of events from 
categories 1 and 4 in that m(S~7r + [ir + ]) range. 

• The fitted yield from the sidebands is subtracted from the fitted yield in the signal region. 

This process suppresses the category 4 background but retains the signal with high efficiency (~ 95% 
at low p* and ~ 90% at high p*). A small fraction of category 4 events have peaking structure in 
m(S~7r + [ir + ]) and survive the sideband subtraction; this rate is 1% or less of the category 1 rate 
in all cases so we neglect it. 

After performing the sideband subtraction described above, the numbers of background-subtracted 
S c + and are 3341 ± 375 and 3195 ± 301, respectively. The p* distributions are shown in Fig. H3 
As discussed in Section Q there are two contributions: S' c from B decays and from the contin- 
uum. S' c produced in B decays have low momentum, especially if the recoiling antibaryon is also 
charmed. Allowing for the motion of the B mesons, the kinematic limit is p* < 2.08 GeV/c, but 
this corresponds to the process B — > S' c p which is heavily suppressed. The limit for Cabibbo- 
allowed processes is p* < 2.02 GeV/c. By contrast, continuum production occurs mainly at higher 
values of p*, with a kinematic limit of p* < 4.63 GeV/c at \J~s = 10.6 GeV. Two separate peaks 
corresponding to these processes are clearly visible in Fig. El 



3.6 Efficiency correction 

For each p* interval, the efficiency e is determined from simulated events in the corresponding p* 
range: 

Yield of true S' c in m(E~7r + [ir + ]) signal window 
Number of generated S' c 

where the yield is obtained by fitting the Am spectrum and true E' c are identified with MC generator 
information. An additional correction is made to take into account signal events which fall into 
the m(S~7v + [n + ]) sidebands and are subtracted from the yield as described in Section 1531 This 
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Figure 6: Background-subtracted S c momentum spectra, not corrected for efficiency. The yield in 
each p* interval is shown for (a) S c + and (b) S c °. 
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Figure 7: Efficiency as a function of p* for (a) S' c + and (b) The factor B(A —> pir ) is not 

included. 

correction, 5e, is obtained from the m(S _ 7r + [-7r + ]) lineshape of simulated S' c in the relevant p* 
range: 

Yield of true E' c in m(S' _ 7r + [7r + ]) sideband 
Yield of true S"^ in m(S' _ 7r + [7r + ]) signal window 

This is a small effect (approximately 1% for E' c produced in B decays and 2% for E' c produced from 
the continuum). The overall efficiency, e(l — r5e), is shown in Fig. [7| 

4 SYSTEMATIC STUDIES 

The following systematic effects are considered. The first four are applied only to the overall 
normalization; the others are treated as fully uncorrelated and are applied to each p* interval 
separately. All uncertainties quoted are relative. 

Particle identification: A systematic uncertainty of 3.5% is assigned to the efficiency of the 
proton identification, as in a previous analysis of the E c system (7j. 
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Tracking efficiency: To correct for a known discrepancy in tracking efficiency between data and 
simulation, systematic corrections of (2.35 ± 7.0)% and (1.55 ± 5.6)% are applied to the 
efficiency for reconstruction of — ► S~tt + tt + and £® — ► S' _ 7r + , respectively. 

Photon efficiency: Based on studies of photon-finding efficiency in control samples, a systematic 
uncertainty of 1.8% is applied to the efficiency. 

A branching fraction: The world-average branching fraction is B(A — > ptt~) = (63.9 ± 0.5)%. 
This results in a 0.8% systematic uncertainty. 

Finite simulation statistics: The statistical uncertainty in the efficiency calculation in each p* 
interval is applied as a systematic uncertainty to the specific data point. This is 5% or lower 
in each interval. 

Signal fitting procedure: The analysis is repeated with a different functional form for the signal 
lineshape, described in Section IA.2I of Appendix EI For each p* interval, the systematic 
uncertainty is taken to be the difference between the efficiency-corrected yields from the two 
functional forms divided by y2. The value depends on the specific p* interval, but is typically 
around 5-10%. 

Background fitting procedure: The background shape is changed from a second-order poly- 
nomial to a fourth-order polynomial and the fit range is increased substantially. For each 
p* interval, the systematic uncertainty is taken to be the difference between the efficiency- 
corrected yields from the two methods divided by \/2. This varies between p* intervals, but 
is typically around 5-10%. In a few intervals with low yields it rises to 30-50%, but is always 
lower than the statistical uncertainty. 

Efficiency correction within a p* interval: If the simulation does not correctly model the p* 
distribution within a p* interval and the efficiency varies significantly across that interval, 
the efficiency may not be predicted correctly. This effect was studied in a previous analysis 
and found to be a few percent or less with the BABAR simulation [Jj. We assign a systematic 
uncertainty of 4% for each p* interval. 

Intermediate resonances in — ► S~tt + tt + : In the simulation, the three-body decay is 
assumed to be entirely non-resonant. However, structure is observed in the Dalitz plot 
distributions in data. We determine the efficiency for two extreme cases: when the non- 
resonant contribution is 0% (e rcs ) and when the non-resonant contribution is 100% (e n onres)- 
The overall efficiency is then changed to (e res + £ n onres)/2 with a systematic uncertainty of 
(|e res — £nonres|)/2. This only affects the S c + mode. The effect is approximately 15% at low 
p* , dropping to zero at high p* as the efficiency becomes more uniform across the Dalitz plot. 

Finite resolution: The reconstructed p* distribution is the true p* distribution convoluted by the 
resolution function. The resolution varies withp*, but is typically 15-20 MeV/c, substantially 
smaller than the bin size (500 MeV/c). We therefore neglect this effect. 

Admixture of cc and BB simulation: The angular distributions of S' c produced in continuum 
events and in B decays differ, resulting in slightly different efficiencies for the two processes. 
This is modelled in the simulation. In general it is unambiguous which process dominates in 
a given p* interval, but in the p* interval 1.5-2.0 GeV/c both may contribute significantly, 
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Figure 8: The efficiency-corrected, background-subtracted p* spectrum for (a) S' c + and (b) S' c °. 
The curve is the simulated continuum distribution described in Section 03 it is fitted to the data 
for 2.0 < p* < 4.5 GeV/c (indicated by the dashed line). 

leading to a slight dependence of the efficiency-corrected yield on the assumed relative pro- 
duction rates. However, the absolute yield in this interval is small and other uncertainties 
dominate, so we neglect this effect. 

After applying these systematic corrections and uncertainties, the p* spectrum shown in Fig. |H] 
is obtained. The inner error bars show the statistical uncertainty (both from data and simulation), 
the middle error bars show the sum in quadrature of the statistical and uncorrelated systematic 
uncertainties, and the outer error bars (where visible) show the sum of all uncertainties in quadra- 
ture. 

5 PHYSICS RESULTS 

5.1 Production rates 

It is clear from Fig. |H] that there is significant E' c production both in B decays and from the cc 
continuum. We separate the contributions of the two processes as follows: 

• B production of E' c for p* > 2.0 GeV/c is assumed to be zero. 

• Continuum production of S' c for p* > 2.0 GeV/c is taken to be the sum of the measured 
yields in each p* interval above 2.0 GeV/c. 

• The data between 2.0 GeV/c and 4.5 GeV/c are then fitted with a suitable function, described 
below. The function is extrapolated down to p* = and the integral over the range 0.0- 
2.0 GeV/c is taken as the continuum production of S' c in that momentum range. 

• The B production of E' c below 2.0 GeV/c is taken to be the sum of the measured yields 
in each p* interval below 2.0 GeV/c less the estimated continuum production in the range 
0.0-2.0 GeV/c 
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Figure 9: Comparison of the efficiency corrected yields in the off-peak sample with those of the 
full data set (off-peak and on-peak combined). The off-peak yields have been scaled up to account 
for the difference in integrated luminosity, and corrected for the small change in the continuum 
cross-section with y^s. Plot (a) shows the S c + and plot (b) shows the S c °. Systematic uncertainties 
and corrections are not included. 

The continuum function is based on the Bowler fragmentation model JH]; tuned to the global 
BABAR data and implemented within the JETSET ^3] generator. Only the amplitude is allowed 
to float in the fit to the data. The fitted function is shown in Fig. |S1 The x 2 /NDF of the fits are 
1.2/4 and 1.8/4 for S c + and S c °, respectively. 

To test the model dependency, we fit the data with a number of other fragmentation models. 
For these crosschecks, we take the parameterizations to be functions of the scaled momentum x p . 
The data above 2.0 GeV/c are well-described by the Peterson model an d by a baryon-specific 
version of the phenomenological model of Kartvelishvili et al. [2111 • The standard deviation of the 
extracted rates when using these three fragmentation models is quoted as the model-dependent 
uncertainty: 

S c + from continuum: 32681 ± 5516 (exp.) ± 4443 (model) 
S'+ from B decay: 77199 ± 7907 (exp.) ± 4443 (model) 
from continuum: 16356 ± 2509 (exp.) ± 1384 (model) 
S' c ° from B decay: 30782 ± 3088 (exp.) ± 1384 (model), 

where the experimental uncertainties combine both statistical and systematic effects. Excluding 
the normalization systematic uncertainties, these correspond to a statistical significance for S' c 
production in B decays in excess of 12<7 for each mode, and a significance for continuum production 
at p* > 2.0 GeV/c in excess of 6a for each mode. As an additional crosscheck, the continuum 
production for p* < 2.0 GeV/c is measured in the off-peak data sample alone. This procedure is 
model-independent (to 2-3%) but has a much larger statistical uncertainty. The results are shown 
in Fig. Inland are consistent with the yields quoted above within statistical uncertainties. 

Dividing the above yields by twice the total number of BB pairs in the data sample, we measure 
the product branching fractions as: 

B(B S'+X) x B(E+ -> = , "vr + vr + ) = [1.69 ± 0.17 (exp.) ± 0.10 (model)] x 10~ 4 
B(B S' C °X) x B{~ Q C -» S-TT+) = [0.67 ± 0.07 (exp.) ± 0.03 (model)] x 10~ 4 . 
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Comparing the second measurement with a previous BABAR result [7], 

B(B -» ~° C X) x £(,~ c -» ~vr + ) = [2.11 ± 0.19 (stat.) ± 0.25 (sys.)] x 10~ 4 , 

we observe that approximately one third of produced in B decays come from E® decays. 

Correcting for the 1/s scaling of the continuum cross-section for data taken at ^/s = 10.54 GeV 
and taking into account the 1.5% systematic uncertainty on the integrated luminosities quoted in 
Section |2 the cross-sections at y/s = 10.58 GeV are: 

a(e + e~ -> ~' c + X) x B(E+ -> = , "vr + vr + ) = 141 ± 24 (exp.) ± 19 (model) fb 
a{e + e~ -> S' C °X) x 5(=' c -» ^-yr + ) = 70 ± 11 (exp.) ± 6 (model) fb. 

Comparing this with a previous BABAR result [7J: 

o-(e + e" -» S°X) x -» S-7T+) = 388 ± 39 (stat.) ± 41 (sys.) fb 

we observe that about 18% of produced in the continuum come from a S'®. 

5.2 Helicity angle distribution 

The quark-model predicts J p = ?j + for S' c , S c and S~ . Under this assumption, the helicity angle 
distribution for the decay processes studied should be flat in the cosine of the helicity angle, cos Oh, 
where Oh defined as is the angle between the S~ direction in the S c rest frame and the S c direction 
in the E' c rest frame. If the E' c has J = |, the angular distribution is of the form c\ + C2 cos 2 Oh 
where c\ and C2 are unknown parameters and C2 may be zero. In general the distribution for higher 
spins is a higher-order polynomial — but the distribution may be fiat if the higher-order coefficients 
are zero. Therefore, a non-fiat distribution in the data would exclude J =\, but a flat distribution 
would not exclude J > |. 

The method used to measure the helicity angle distribution in the data is similar to the one 
used for the p* spectrum: we divide the data into six slices of cos Oh and fit the mass spectrum in 
each slice. A p* threshold of 2.5 GeV/c is applied throughout to improve the signal-to-background 
ratio. After correcting for efficiency, the distributions are fitted with two functions: first, a flat 
distribution: 

fo(cos h ) = a (3) 

and second, a symmetric quadratic: 

/o(cos^) = a (l + (3 cos 2 Oh) . (4) 

Separate fits to each mode in the data are performed and the fit results are shown in Table ^ 
The data are clearly consistent with being flat (x 2 /NDF less than unity). The fitted quadratic 
parameters (3 are consistent with zero, though with large statistical uncertainties. Since the two 
modes should have identical helicity distributions, we weight them according to their statistical 
precision and combine them — this is shown in Fig. ^] and the fit results are given in Tabled The 
data are still fit well by a flat distribution, though with a somewhat reduced x 2 probability (20%). 
From this we conclude that the data are consistent with the predicted J = \ but that higher spins 
cannot be excluded. 
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Figure 10: Fits to efficiency-corrected helicity distributions for data with p* > 2.5 GeV/c. The 
plot shows the normalized, weighted sum of the S c + and S c ° distributions. The solid line assumes 
a flat helicity distribution, whereas for the dashed line a quadratic term is added. 



Table 1: Results of fits to efficiency-corrected helicity distributions. The x 2 goodness-of-fit and 
the number of degrees of freedom are given for the flat and quadratic distributions given in Eq. |3] 
and along with the fitted parameter j3 from the quadratic distribution. The results are given for 
the S c + and S" c ° samples individually, and for a weighted sum of the two samples. 

S' c + S' c ° Weighted Sum 

X 2 /NDF for flat 4.0/5 (55%) 4.4/5 (50%) 7.3/5 (20%) 

X 2 /NDF for quadratic 2.7/4 (61%) 1.3/4 (87%) 3.0/4 (55%) 

for quadratic 0.63 ± 0.68 1.04 ± 0.82 0.79 ± 0.52 
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6 SUMMARY 



We have confirmed the CLEO observation of the S' c + and S'® states from the cc continuum. In 
addition, we found that B mesons decay at a substantial rate to S c + and S c °. This is the first 
observation of such decays; the statistical significance is in excess of 12a for each mode. We have 
measured the production rates of S' c from B decays (expressed as a branching fraction) and from 
the cc continuum (expressed as a cross-section); in both cases, the absolute rate is scaled by the 
unknown absolute S c branching fraction. We have measured the angular distribution of S' c — > S c j 
decays and found it to be consistent with the prediction for J p = | + . However, higher spins cannot 
be ruled out. 
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Appendix 



A Lineshape parameterizations 

A.l Triple Gaussian function 

The lineshape function /(Am) is parameterized as follows: 



/(Am) = N(l - h ~ /s)G(Am; Ml, a x ) + Nf 2 G(Am; M i + A/i 2 , <r 2 ) + Nf 3 G(Am; M i + A/i 3 , a 3 ) (5) 



where G(x; /x, cr) is a Gaussian function of unit area with mean [i and width a. The nine parameters 
are then interpreted as follows: 



N Total fitted yield 

fi Fraction of yield in second Gaussian function 

/ 3 Fraction of yield in third Gaussian function 

fi\ Signal mass 

Afi2 Mean of second Gaussian function with respect to the signal mass 

A^3 Mean of third Gaussian function with respect to the signal mass 

o\ Width of first Gaussian function 

cr 2 Width of second Gaussian function 

cr 3 Width of third Gaussian function 



When fitting simulated events to determine the lineshape, all nine parameters are allowed to vary 
independently. In order to improve fit convergence, the following bounds are placed on the variation 
of the parameters: 



N No bounds 

fi No bounds 

/3 No bounds 

Hi Limited to (—5, +5) MeV/c 2 relative to the true mass 

A/i2 Limited to (—15, +5) MeV/c 2 relative to \i\ 

A/i3 Limited to (—35,-5) MeV/c 2 relative to /xi 

ax Limited to (0, 10) MeV/c 2 

cr 2 Limited to (4,20) MeV/c 2 

o- 3 Limited to (20, 100) MeV/c 2 



A. 2 Alternative lineshape parameterization 

The following functional form was also used as a cross-check and to determine the systematic 
uncertainty due to the signal lineshape: 
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where 



■r-2 



Xp+ 2 



£_ 



+ 1 



Zip controls the 



Of the six parameters A p controls the amplitude, x p controls the peak position, 
width, £ controls the asymmetry in the central region, p\ controls the lower tail, and P2 controls 
the upper tail. 



B Individual mass spectra 

The fitted mass spectra for individual p* intervals are shown in Fig. ^2andElf or " C + and E®, 
respectively. 
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Figure 11: The £^7 invariant mass spectra for each of the ten p* intervals used. The points show 
data from the signal mass region (within 3<r of the central value). 
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Figure 12: The ^7 invariant mass spectra for each of the ten p* intervals used. The points show 
data from the signal mass region (within 3cr of the central value). 
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